Intra- site 4/ — 5d electronic correlations in the quadrupolar model of the 7-a phase 

transition in Ce. 
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As a possible mechanism of the 7 — a phase transition in pristine cerium a change of the electronic 
density from a disordered state with symmetry Fm3m to an ordered state Pa3 has been proposed. 
Here we include on-site and inter- site electron correlations involving one localized 4/-electron and 
one conduction 5d-electron per atom. The model is used to calculate the crystal field of 7-Ce and 
the temperature evolution of the mean-field of a-Ce. The formalism can be applied to crystals where 
quadrupolar ordering involves several electrons on the same site. 
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I. INTRODUCTION 

Elemental solid cerium is known to undergo struc- 
tural phase transitions BB In the pressure-temperature 
phase diagram of Ce the puzzling long-standing problem 
is to understand the apparently isosteuctural transition 
between the cubic 7 and a phases .lilt! An isostructural 
phase transformation can not be ascribed to a conden- 
sation of an order parameter and therefore can not be 
explained by the Landau theory of phase transitions. In 
the past, several models and-thcorics have been suggested 
to address this problemJjEj. Among them, a Mott like 
transition, fee 4/ electronsoiJ and a Kondo effect based 
approachErtHI are the competing ones. Also, new compu- 
tational schemes have been applied to the problem us- 
ing dynamical mean-field, theory combined with the local 
density approximation£3 

Under pressure above 5 GPa a-Ce becomes unstable 
and transforms first to a crystal with C2 jm or a-U space 
symmetrya (a'-Ce) and then to a body-centexed tetrag- 
onal (bet) structure (a"-Ce) above 12 GPa.EJ This se- 
ries of transformations can not be explained by invoking 
the concepts of the 4/ localization-delocalization tran- 
sition or Kondo volume collapse models and indicates 
that there are anisotropic interactions present in the a 
phase of cerium. Such electron interactions can be of 
quadrupolar origin that are known to drive symmetry 
lowering phase transitions in many lanthanide and ac- 
tinide compounds .113 

Recently, the isostructural character of the 7 — a 
phase transition has been questioned by Eliashberg and 
CapellmannEa They suggest that a-Ce should have a dis- 
torted fee structure. Independently, the present authors 
have p ut ., fo rward a theory of quadrupolar ordering in 
cerium.EZHi3 There, it was proposed that the 7 — a trans- 
formation is not really isostructural. Rather, it was-.as- 
sociated with hidden electronic degrees of freedomO In 
our previous work, Refs. 17 Tq, we have suggested that 
the symmetry change is from Fm3m to Pa3. This sym- 
metry lowering is a special one. Although accompanied 
by a lattice contraction, it conserves the fee structure of 



the atomic center of mass points (cerium nuclei) and is 
solely due to orientational order of electronic densities. 
Such a scenario reconciles the 7 — a transformation with 
the Landau theory of phase transitions. Our consider- 
ations for cerium have been inspired by the theoryEEl of 
orientational phase transition in solid Cgo where a sim- 
ilar space symmetry change (FmSm — > Pa3) occurs at 
255 K at room pressure.Ell 

The present article is a continuation of our approach 
to the problem of the 7 — a transition in Ce based on 
the technique of multipolar interactions between elec- 
tronic densi ties of conduction and localized electrons in a 
crystal J 1 Tr 8 ! Our second motivation is to extend our ini- 
tial method, Ref. [l?], for the case when two electrons 
(/— and d— ) are at the same site of cerium and all 
intra- site interactions (including the on-site exchange) 
between them are taken into account. Our treatment of 
intra- site correlations is closely related with the method 
used bv-Condon and Shortley for many electron states of 
atoms.E3 Provided that the average number of electrons 
per site is conserved this method is exact for intra- site 
correlations and -goes beyond the usual self-consistent- 
field approaches c3 employed by band structure calcula- 
tions. 

Besides the problem of the 7 — a phase transition m 
Ce, the microscopic method can be applied furtherE£l 
to describe quadrupolar ordering and to perform crystal 
field calculations of many / electrons on the same site. 
There are numerous compounds exhibiting quadrupolar 
ordering at low temperatures^ and there is a sustained 
interest in undecstandiiig. their properties. |-Thus, Efe 
centlv-pyB 2 C 2 ,tZl DyB 6 Ji UCu 2 SnJ^ PrPb 3 ,El YbAs,H 
YbSbE3 were reported to undergo a quadrupolar order- 
ing. 



II. THE MODEL 

As follows from electronic band structure calculations 
of 7-Ce there exist three conduction electrons per atom 
which form the (6s6pbd) 3 metallic band and one local- 
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ized 4/ electron.ElEH] In our previous work, Ref. [l8[ we 
have already considered electric multipole interactions 
between conduction electrons and the localized 4/ elec- 
trons. Below we focus on the on- site and inter- site 
correlations in the system and will simplify the model. 
We consider the instantaneous configuration 6s 2 5g?4/ as 
having the largest statistical weight on a cerium site (in 
comparison with other possibilities such as 6s6p5<i4/, 
6s 2 5d 2 , 6s 2 4/ 2 and etc.). The two 6s electrons give only 
a spherically symmetrical density on the cerium site and 
their lowest energy state corresponds to a singlet. There- 
fore, as a first approximation we discard them as giving 
a closed shell. (Indeed, the level structure of 6s 2 5c?4/ 
corresponding to atomic Ce I is similar to that of 5d4/ 
of La IIo) In fact, in doing so we omit the s — d elec- 
tron transitions which can contribute to the quadrupolar 
density£3 We are left then with one 5d conduction elec- 
tron and one localized 4/ electron. In the electronic band 
structure calculations the charge density of the 5d elec- 
tron on the cerium site is considered as an average over 
all occupied k states (E(k, a) < Ep, where a is the band 
index and Ep is the Fermi level). We have shown in Ap- 
pendix B of Ref. |l7| that the electron density is mainly 
spherical which corresponds to the standard "muffin-tin" 
(MT) treatment in electronic band structure calculations. 
The spherical density of 5d and 4/ electrons will be the 
starting point in this work. We consider the 5d electron 
on a cerium center being instantaneously coupled with 
the 4/ electron and include in the model all correspond- 
ing intra- site interactions, crystal electric field effects 
and inter- site multipolar electric interactions. From the 
technical point of view, this /d-model is a many electron 
generalization of the concepts of Ref. [l7|. We are aware 
that the model based on the /d-configuration is incom- 
plete, but it has an advantage of taking into account all 
intra- site interactions (often referred to as Hund's rules) 
which are usually omitted in the electron band structure 
calculations .Bo Later we will briefly discuss a possibil- 
ity to refine our model with the help of the valerppJbond 
(or Heitler-London) theory of chemical bondingH'Ej 

We consider a face centered cubic (fee) crystal of N 
Ce atoms. Each atomic site possesses one 4/ and one 
5d electron. The position vector of an electron near a 
crystal lattice site n is given by 



R(n) =X(n) + r(n). 



(2.1) 



Here X(n) is the lattice vector which specifies the cen- 
ters of the atoms (or Ce-nuclei) on a rigid fee lattice. 
The radius vector r(n) is given in polar coordinates by 
(r(n), fi(n)), where r is the length and O = (6, 4>) stands 
for the polar angles. We label the two-electron basis ket- 
vectors at a lattice site n by a single index Ifd or, alter- 
natively, by the pair of single electron indices (if, i d )- 



\I 



fd)n 



(2.2) 



The index i stands for the electron orbital and spin pro- 
jection quantum numbers. The corresponding basis wave 



functions are 

(r,r"\Ifd)n = (r\if)fi ' (r'\i d )n, 

where 

(f\i f ) H = K f (r(n))(h\i f ), 
(f'\i d )rt=Tl d (r'(fi))(n'\i d ). 



(2.3) 

(2.4a) 
(2.4b) 



Here IZf and TZ d are radial components of the 4/ and 
the 5d electron, respectively; h stands for Cl(n). There 
are 14 orientational vectors (or spin- orbit als) (n\if) for a 
4/ electron (z y =1-14) and 10 orientational vectors (n\i d ) 
for a 5d electron (^=1-10). These spin-orbitals can be 
written as 



(h\i f ) = (n\m f }u s (s z (f)), 
(h\i d ) = (n\m d }u s (s z (d)). 



(2.5a) 
(2.5b) 



Here u s is the spin function (s = ±) for the spin projec- 
tions s z — ±1/2 on the z-axis. The orbital parts, (n|m/) 
(m/=l-7) and (h\m d ) (m<i=l-5), are expressed in terms 
of spherical harmonics YJ m (f2) = (n\l,m). We find it 
convenient to work with real spherical harmonics. We 
consider 



{y 3 °, r 3 lc , y 3 1s , y 3 



2c 



, Y 3 2s , Yt, Y 3 3s } = (h\m f ) (2.6a) 



for 4/ electron (corresponding to m/ = 1 — 7) and 



\rlc 
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Y 2 ls , Y 2c , Y? s ) = (h\m d ) 



(2.6b) 



for 5d electron (corresponding to = 1 — 5). We use the 
definition of real spherical harmonics of Ref. [37] (see also 
explicit expressions (2.1) in Ref. |l7|) that-is different from 
the definition of Condon and Shortley.E3 The advantage 
of using the basis with real spherical harmonics is that the 
matrix elements of Coulomb and exchange interactions 
stay real. 

The order of indices in @.2| ) and ( |2.3| ) is important if 
we associate the first electron with the / state if while 
the second with the d state i d . Then in addition to the 



vectors (2.2) we have to consider the states described by 
the vectors \i d ; if)n (the first electron is in the i d state 
and the second is in the if state). However, from the dy- 
namical equivalence of the electrons we can permute the 
spin-orbitals to the standard order, Eq. (p.2|), by using 



lid! lf)n 



(2.7) 



since it requires the interchange of the two electrons. In 
order to describe the sam e quantum st ate (if, i d ) we will 
use the basis vectors (2.2) and apply (|]^) when needed. 
(Alternatively, one can use the procedure of antisym- 
metrization of the basis v ecto rs ( [2.2] ) as described else- 
where.) Thus, our basis (2.2) consists of 140 different 
vectors \If d ). 

The ground state energy of the (4/5d) electron sys- 
tem, Eq, can be calculated in local density approxima- 
tion (LDA) with spherically symmetric Coulomb and ex- 
change potentials. Going beyond this model in atomic 
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cerium, one has to take into account multipolar on-site 
(also called intra- site) Coulomb interactions and spin- 
orbit coupling. In solid cerium, the interactions with 
conduction electrons and inter- site Coulomb interactions 
still have to be added. 

In the following we will study these effects within 
a unified formalism based on a multipole expansion of 
the Coulomb potential and of the systematic use of site 
symmetry of the crystal lattice. For the case of on- 
site Coulomb interactions between two electrons (charge 
e = —1) we have 



V(R(n),R'{n)) 



1 



f'(n)\ 



(2.8) 



The multipole expansion in terms of site symmetry 
adapted functions (SAF's) S A (n) reads: 

V(R(n),R'(n')) = £ v AA (r, r') S A (h) S A (n% (2.9a) 



where 



v AA >{r, r') 



4tt 



-('+!) 21 + 1 



>AA' 



(2.9b) 



with r > = max(r,r'), r< 



min(r,r') and i5aa' 



S tt i5u'. Clearly, the last expression is site independent. 
The SAF's are linear combinations of spherical harmon- 
ics and transform as irreducible representations of the 
site point group, Ref. [?7|. The index A stands for (Z,r), 
with t = (r, /i, k). Here I accounts for the angular de- 
pendence of the multipolar expansion, T denotes an irre- 
ducible representation (in the present case of the group 
Oh), [i labels the representations that occur more than 
once and k denotes the rows of a given representation. 

On the other hand, the Coulomb interaction between 
two electrons at different sites n ^ n' (inter- site) reads 



1 



V(R(n),R'(n r )) = — 

\R(H)-R'{n')\ 

The multipole expansion is given by 



(2.10) 



V(R(n),R'{n')) = V>aa' r, r') S A {n) S A > (ft'), 



AA' 



where 



v AA >(fi,n'; r,r') = 



dfl(n) J dfi'(n') 



(2.11a) 



S A (n) S A >(n') 
\R(n)- R!(n')\ 

(2.11b) 



The inter- site multipole-expansion (2.11a) is anisotropic 
and converges fast sinceE3 



v AA ,(ri,ri'; r,r') 



(r)V) 



\X(n) - X(n')\ 



(2.12) 



Therefore, it is sufficient to consider it only for nearest 
neighbors. From the practical point of view, one can 
calculate v AA i (ft, ft'; ro, r' ) only for two fixed radii ro and 
Then one obtains v AA i (ft, n'; r, r') as a function of r 



and r' by employing the dependence ( 2.12| ) 



III. THE INTRA- SITE INTERACTIONS 



The interactions which we analyz 
present already in atomic cerium.' 22 ' 



n this section are 
We have consid- 



ered a part of these correlations in previous wor kjll and 
here we study all of them in the framework of the fd— 
model. In fact, these multipole interactions are responsi- 
ble for the electronic terms of atoms£3 We will see that 
their combined effect lowers the energy of the cerium 
atom by ~ 1 — 2 eV in comparison with the spher- 
ically symmetric case, yet usually they are not taken 
into account in the electronic band structure calcula- 
tions in solids. Although our consideration in this sec- 
tion is based-, on the original technique for multipole 
interactions, Ell it overlaps largely with the method of 
Condon and Shortley.EHl However, for two reasons we have 
decided to briefly review it here. Firstly, we consider be- 
low a more general case which is not limited by the LS 
(Russell-Saunders) coupling and the consideration of di- 
agonal matrix elements. Secondly, the results of this sec- 
tion are used to describe crystal electric field effects and 
the phase transition to the Pa3 structure. 

The direct matrix elements for the intra- site Coulomb 
interactions are obtained if we consider only the / — / 
transitions for the first electron and the d — d transitions 



for the second. We start from Eq. (2.9a) and obtain 



(lMnV(R(n),R'(fi))\J fd )%° ul 
= U AA c A{ifif) c A {i d jd), 



where 



; FD 

J AA 



(3.1) 



drr 2 dr 1 r"n}{r)n 2 d {r')v AA {r,r') (3.2) 



accounts for the avera ge r adial dependence while 
VAA{ r , r ') is given by Eq. (2.9b). We use the superscripts 



F and D in order to indicate that we have transitions 
between two 4/ states [F = (/,/)) and the transitions 
between two 5d states (D = (d, d)). The elements c A are 
defined by 



(3.3a) 



c\{ifjf) = / dCl(i f \h)S A (n)(h\jf), 



c A {idjd) = I dVt (id\n)S A (h)(h\jd). 



(3.3b) 



The other possibility is to consider the transitions 
4/ — > 5c? for the first electron and the transitions 5d — * 4/ 
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for the second. This give s the exchange interactions and 
then we should use (2/7) in order to return to the stan- 
dard order of the spin-orbitals. We find 

(I fd \ ri V(R(n),R'(n))\J fd )r h 

= -E U AA )W)c A(vi d )cA(W/), (3.4) 
A 

where 

v (fd)(df) = J drr 2 j ' dr , r > 2 n f (r)K d (r)K d (r')Kf(r') 
xu AA (r,r'), (3.5) 

and 

CA(ifjd)= JdQ{if\h}S\{n){h\jd), (3.6a) 



(idjf)= dn(i d \h)S A (n)(n\j f } 



(3.6b) 



We o bserve that in the basis with real orbitals, Eq. 
(2.6a,b), and with real functions S\ the coefficients ca 



are real and we get 



c\(idjf) = c\{jfi d ). 



(3.7) 



We start with the description of spherically symmetric 
terms (/ = 0) corresponding to the trivial function Sq = 



1/V47T. The coefficients ca in (3.3a,b) become diagonal 



ci=o(idjd) 



1 



47T 



c i=o(if3f) 



(3.8) 



while Ci—o(itjd) = ci=o(idjf) = 0- Hence, we obtain a 
contribution to (I\V\J} Coul which is proportional to the 
unit matrix. Since it corresponds only to a shift of the 
ground state energy, it is irrelevant. 

In consi derin g the other contributions (with I > 0) in 
(3.1) and (3.4) we will take advantage of the selection 
rules imposed by the coefficients ca, Eqs. ( |3.3a| ,b) and 
( 3.6a ,b). First of all, we notice that the coefficients ca 
are diagonal in terms of spin components u s . From the 
theory of addition of angular momenta (see, for example, 
Rcf. |36| ) we know that nonzero coefficien ts ca (i fj f ) can 
occur if I (in A) equals to 0,1, 2,.. .,6, Eq. (3.3a). Further- 
more, the odd va lues o f I are excluded due to the parity of 
the integrand in (3.3a) and finally we obtain that I — 0, 2, 
4 and 6. Analogously, for the d—d transitions the allowed 
coefficients are with I = 0, 2 and 4, Eq. ( (3.3b| ). For the 
f — d transitions, Eq. ( |3.6a ,b) we find that I = 1, 3 and 
5. Next we notice that if the radial parts, IZf, 7Z d , are 
the same for all spin-orbitals of 4/ and 5d sta tes, cor- 
respondingly, then the integrals (3.2) and (3.5) depend 
only on I. We can condense the notation, 



(fd)(df) 
A A 



AA 



for 



for I' = 1, 3 and 5. 



I = 0, 2 and 4; and v\I ' = v 
In fact, these integrals are proportional-to F k and G k 
the notations of Condon and Shortley,c3 



.FD 4?r z-2 



.Xfd) _ 



—F\±f,bd), 

47T 



.FD 47r it ! 



^G 1 (4/,5d), vi fd) = ^G 3 (4/,5d), (3.9) 



Y^ 4 (4/,5d) 

47T 



.Xfd) 



4ir 



= — G 5 (4/,5d). 



(Notice however that our coefficients ca are different from 
those of Co ndon and Shortley.t3) We have estimated the 
integrals in ( |3.9| ) from the radial dependences IZf and lZ d 
obtained from calculations on a cerium atom in local den- 
sity approximation (LDA). We have found that wf D = 

85858, < D = 24260, v[ fd) = 83128, v^ fd) = 26563 and 



Xfd) 



12584, in K (Kelvin). 



Finally, we rewrite expressions (3.1) and (3.4) in ma- 
trix form as 



(If d \V{intra)\J fd 



Coul 



FD FD 



(I fd \J fd ) + v[ D c! D (I fd \J fd ) 



(3.10a) 



and 



(Ifd\V(intra)\Jf d 



, exch 



[v[ fd) c{ fd \lfd\Jfd) 



+v 



(fd)(fd) 



(Ifd\Jfd)+vi fd) ci fd) (I fd \Jfd)}, (3.10b) 



Here the direct Coulomb matrices cf D (If d \Jf d ) are de- 
fined as 



„fd 



(Ifd\Jfd) = (^C(!,T)(if3f)c(i,T)(idd^j , (3.11a) 



where I = 2 and 4; and the three "exchange" matrices 
c j \ (Ifd\Jfd) (Z = 1, 3 and 5) are given by 

C^HtfdlJfd) = ^C{l,r){ifjd)C{l,r){idif) ■ (3.11b) 



We solve the secular problem for the 140 x 140 matrix of 
intra- site interactions 



(I f d\V(intra)\J f d) = 
(I fd \V( ln tra)\J fd ) Coul 



{I fd \V(intra)\J fd ) exc \ (3.12) 



and obtain the ten energy levels Ef d quoted in column 2 
and 3 of Table I. This term spectrum corresponds to the 
usual LS (Russell-Saunders) coupling without spin-orbit 
interactions. We have als o che cked our results by working 
with the eigenvectors of (3.12) which can be obtained in- 
dependently by exploiting the formulas (Table 4 3 of Rcf. 
|22] ) of vector addition of angular momenta (j\ — 3 and 
j2 = 2). The levels of Table I correspond to the followiug 
parameters in the notations of Condon and ShortleyEHl 
F 2 = 325 A, F 4 = 25.1, Gi = 567, G 3 = 47 and G 5 = 7.2 
(in K) (These parameters should not be confused with 
those in (p|)). 

We now consider the effect of the spin-orbit coupling. 
Starting with the spherically symmetric LDA calculation 
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TABLE I. Term energies of the fd configuration in the 
absence of spin-orbit coupling, E m and AE stand for the sin- 
glet-triplet energy means and the singlet-triplet energy differ- 
ences, respectively. L is the orbital quantum number of the 
two electrons. The energy corresponding to the spherically 
symmetric description is taken as zero. All energies are in 
units K. 



L 


singlets 


triplets 




AE 


P 


13541.8 


5384.5 


9463.1 


8157.3 


D 


-3011.8 


1951.8 


-530.0 


-4963.5 


F 


2767.8 


-6616.1 


-1924.1 


9383.9 


G 


-9378.4 


-1485.3 


-5431.8 


-7893.1 


H 


12310.8 


-5653.4 


3328.7 


17964.2 



TABLE II. Calculated lowest term energies of the fd 
configuration with the spin-orbit coupling and the Lande 
<?-factors (columns 2-4), AE — Efd — E{d( l Gi). Last two 
columns are experimental data, Ref. 34. All energies are in 
units K. 





Efd 


9 


AE 


Ce I 


La II (fd) 


X G 4 


-10570.3 


0.9323 











3 F 2 


-9226.5 


0.7030 


1343.8 


329.2 


881.7 


3 H, 


-8155.1 


0.8958 


2415.2 


1840.9 


1764.7 


3 F 3 


-7048.8 


1.0825 


3521.5 


2393.0 


2354.5 


3 H 5 


-6449.4 


1.0334 


4120.9 


3177.9 


2850.7 


3 G 3 


-4247.3 


0.7659 


6323.0 


1998.5 


5472.8 



of a cerium atom we obtain that in the one-electron ap- 
proximation A so (/) = E f (7/2) - E f (5/2) = 4003.4 K 
and A so (d) = E d (5/2) - E d (3/2) = 2344.0 K. This gives 
for the spin-orbit coupling constants Q = 1143.8 K and 
Q = 937.6 K. Therefore, a typical value of spin-orbit 
splitting is ~ 1000 K which shows that it can not be 
treated as a small perturbation to the LS term scheme, 
Table I. In order to take into account the spin-orbit cou- 
pling exactly we have to consider the operator 



V so (f) + V so (d), 



where 



V so (f) = ( f L(f)-S(f), 
V so (d) = QL(d)-g(d). 

The matrix of interactions reads 



(3.13) 

(3.14a) 
(3.14b) 



(If d \V so \J fd ) = {(if\V so (f)\j f ) 6 idjd + (i d \V so (d)\j d ) S ifjf ) . 

(3.15) 

Since we know the matrix elements (if\V so (f)\jf) and 
{id\V so {d)\jd) (see for example the explicit Eq. (A. 2) of 
Ref. |l|) we can calculate the matrix elements for the fd- 
configuration. We now solve the secular problem for the 
sum of the intra-site and spin-orbit interactions, starting 
from the matrix of 



U\ mtra = V{intra) + V SOl 



(3.16) 



and obtain 20 energy levels {Efd}. Since we are inter- 
ested only in the lowest levels, we quote in Table II the 
first six out of the 20 levels. (The procedure of calcu- 
lation of magnetic moments and the Lande g-factors is 
given in Appendix A.) ._. 

Comparison with the experimental datao on a Ce 
atom shows that the order of the three lowest levels, 1 G < 4, 
3 F 2 and 3 i?4, with the multiplet 1 G 4 as the ground state 
is correct. On the other hand, the experimental data 
show that in atomic Ce different levels (such as 1 G, 3 F, 
3 H) are considerably mixed up. Due to the strong spin- 
orbit coupling there is an appreciable mixing between 



1 G , 4 and 3 i/ 4 , and between 3 F 2 and 1 D 2 . The ground, 
state of a Ce atom has 55% of X G± and 29% of 3 H 4 S3 
The next level which lies only 329 K above the ground 
state has 66% of 3 F 2 and 24% of 1 D 2 S3 In conclusion, 
although the actual atomic spectra of cerium differ some- 
what from our results obtained for the fd— configuration, 
Table II, our approach captures the main properties and 
we will use it for Ce in the solid state. In the following we 
will extend the calculations of the energy level scheme by 
including inter-site interactions. We will treat separately 
the 7 and the a phase of solid Ce. 



IV. INTER- SITE INTERACTIONS 

In this section we will first consider the matrix of inter- 
site interactions for the (fd) system on a crystal in gen- 
eral. Next we will show that the interaction is largely 
simplified by crystal symmetry and derive the crystal 
electric field (CEF). We calculate the energy spectrum 
of the (fd) system in presence of the crystal electric field 
in the disordered cubic phase. 

We start from expression ( |2.11a ) and write it in the 
space of two electron state vectors \If d )- Carrying out 
the angular integrations dVl(n), dO'(n), dtt(ri r ), dQ'(ri r ), 
we obtain 



a/3 a'P' AA' 

X {cA(iaja) S(i/3j/3)} {CA'(i a >ja') 8(ip>jp>)} , 



(4.1) 



Here each of the indices a, a', (3, (3' runs over the l abels 
f, d. The coefficients c\ are defined by Eqs. ( 3.6a ) and 
( |3.6b| ), while ^(ipjp) stands for the Kronecker delta sym- 
bol. For a = /, we have [3 = d and for a = d, f3 = f, 
with a similar correspondence between a' and (3'. The 
inter- site interaction element is given by 

vl a % a ' (n - H') = J dr r 2 J dr' r' 2 

X 1Zl(r)K 2 a ,(r')v AA ,(n,n'; r,r'), (4.2) 
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Notice that only direct Coulomb interactions are 
present in ( |4.l| ) , and with the help of the selection rules 
for the coefficients c\ we conclude that only the interac- 
tions with I = 0, 2, 4 and 6 have to be considered. 



A. The crystal electric field (7-phase) 

In the 7-phase the electronic density is compatible with 
the crystal structure Fm3m. At each atomic site the 
crystal electric field (CEF) has the point group symme- 
try Oh- In lowest approximation, the CEF corresponds 
to the potential experienced by a charge at a central 
site n, when spherically symmetric {V — 0) contributions 
from charge densities at the twelve neighboring sites n' 
on the fee lattice and similar terms from the electronic 
density in thpJnterstitial regions are taken into account. 
PreviouslyOES these effects were studied for a single 4/ 
electron per Ce atom. Here we present an extension to 
the (4/5d) system. 

In crystal field approximation the functions SA'(n') at 
any of the twelve sites ft' reduce to Y§ = 1/ v4~7r". We will 
write an index for A' (I 1 — 0, A\ g ). The coefficients ca' 
in Eq. (4.1) now reduce to 



Co(ia'ja> 



(4.3) 



At the central site n, the electronic density has full cu- 
bic symmetry. We denote the corresponding SAF's by 
Sa.1 (n), Ai = (/, A\ g ), where A\ g stands for the unit rep- 
resentation of the cubic site group Oh- We retain the 
functions for I = 4 and 1 — 6, which correspond to the 
cubic harmonics K4 and K e . The selection rules imply 
that the d — d transitions are perturbed by K4 only, while 
for the / — / tran sitions both K4 and _Kg are relevant. 
Expression (2.11a) reduces to 



V(R(H),R'(n')) = --Lj^ofafJ'; r, /) S Al (n). 

(4.4a) 



Ai 



The elements 

v Al0 (ri,ri'; r,r') 
1 



\R(n) - R'(n')\ 



(4.4b) 



have the same value for all twelve neighbors ft' on the fee 
lattice. In additio n the y are independent of r' , as follows 
from expression (2.12) for V = 0. We then define the 
crystal field operator by 



V CF (R(n)) 
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y^i> Al o(n,n'; r, /') S Al (r 



(4.5) 



Ai 



Returning to expression (4.2) we write within the crystal 
field approximation 



vt% a (H-n')=v^' -Q a/ , 



where 



v li'= / drr n a ( r ) v Mo(n,fi; r, / ) 



and 



dr' r ,2 TZ 2 a ,(r r ). 



(4.6a) 



(4.6b) 



(4.6c) 



Notice that Q a > stands for Qf or Qd which are 
(in units e) of the 4/ or 5d electron. As before 
integration is taken over < r 1 < Rmt, where Rmt is 
the radius of the muffin-tin sphere. Besides 4/ and 5d 
electrons we can also consider similar contributions from 
6s electrons and nuclei belonging to nearest neighbors. 
Noti ce th at the interaction parameters vj^Q and Vj^', 
Eq. flOg ), remain the same for all these contributions 
and all we have to do is to collect the charges together. 
Finally, after summation over 12 nearest neighbors and 
simplifications, we obtain 

(Ifd\nVcF(R(n))\J fd ) n = 

[ B L c Ai (ifjf)S(idjd) + B d Ki c Al {i d jd) S(i f 3f)] , (4.7) 

Ai 

where 



B{ 



si 



-j^Qeffev^ 
12 



-in 



(4.8a) 
(4.8b) 



Here again we write D for (dd) and F for (//). We take 
as an effective charge Q e ff = Qmt, the total charge in- 
side p^MT-sphere. In contradistinction to our previous 
workEOlla here we have neglected the effect of intersti- 
tial charges. Previously it was found that a homogeneous 
distribution of negative charge in the interstices increases 
the effective charge by an amount 2.85Qa/t- This fact 
due to the angular dependence of the leading cubic har- 
monic Sad Ai = (I = 4,Ai g ) in Eq. ( |4.4b[ ). Indeed 
K±(fi) is positive and maximum along the cubic direc- 
tion [100] (the centers of the interstices) and negative 
and small along [110] (the sites n'). On the other hand, if 
we consider an inhomogeneous charge distribution where 
most of the electronic density in the interstices is located 
close to [110], then the contribution to the crystal field 
from interstitial charges can be assumed to be negligi- 
bly small. We observe that previously the inclusion of 
contributions from a homogeneous charge distribution in 
the interstices has led to an overestimation of calculated 
crystal field, splitting in comparison with the experimen- 
tal values Jia 



In practice, it is convenient to calculate v^ * f ror 



y AxO 



VAiojn, ft'; R mt, /') a 



R 



(4.9a) 



MT 
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where 



Qi 



dr''. 



■ C+2) 



Kir') 



(4.9b) 



Then the CEF operator (4.5), with an effective charge 
Q e f f , can be rewritten as 



VbiK£(n)) = ^B Al S Al (r 



(4.10a) 



whereE3 

-Bai 
Our 



12 o(H, ft'; Rmt, /') 

Q e // e — 



47T 



calculations 



i? 



for 



,7-Ce 



(4.10b) 

(a=9.753 
=0.22271, 

l\ IT = 0.01739 (ina.u.), and 
fl| =1198.4, fi{ =193.9, 5^=74.4 (all in K).@ We then 
consider the Hamiltonian 



a.u.) yield Qmt = +0.9136 |e|,§ qf/R 4 



MT 



H-<(ri)=U\ mtra + V CF (n), 



(4.11) 



which we associate with the 7 phase of Ce. By diagonal- 
izing H 7 we have found that in the cubic CEF the twenty 
atomic- like levels of cerium are split into 58 distinct sub- 
levels which can be labeled by single valued irreducible 
representations A^T^, E{T 3 ), 7i(r 4 ) and T 2 (T 5 ) oiO h . 
In particular, three lowest levels-of-jcerium are split ac- 
cording to the following scheme,p3'E!3 



: G 4 ^ A l +T 1 + E + T 2 , 
T 2 + E, 

T 2 + A 1 +T 1 + E. 



3 F 2 



3 H 4 



(4.12a) 
(4.12b) 
(4.12c) 



The calculated splittings of these levels are quoted in 
Table III. In presence of a magnetic field, there occurs 
an additional splitting of the triplets. The correspond- 
ing magnetic moments are given in the last column of 
Table III. (Details of the calculations can be found in 
Appendix A.) In the second part of the present section 
we will study the energy levels in the ordered phase. 



B. Quadrupolar ordering (a-phase) 

In the cubic phase (7-Ce) with the Fm3m space sym- 
metry the nontrivial electron density distribution is given 
by cubic harmonics Ki(fl) and K§(£1) with I = 4 and 
6. All quadrupole densities with / = /' = 2 average to 
zero. Only fluctuations of electric quadrupoles are al- 
lowed in the interaction Eq. (4.1) and those lead to an 
effective attractive interaction at the X point of the Bril- 
louin zone.t2ML3 This interaction drives a transition to a 
new phase which is characterized by an ordering of elec- 
tric quadrupoles such that the space group symmetry 
is Pa3. This order-disorder transition is accompanied 



TABLE III. Lowest levels of the energy spectrum of 
H 1 = U\intra + Vc f , J-Ce. Numbers in parentheses stand 
for degeneracy; Aei =1409.6 K, Ae 2 =2499.4 K; the site 
group is O h . 



r // 
1 , n 




ti, 111 IV 


\&i tij, m iv 


Jvlz, 111 flB 


4, i 
j±i, 1 


( L ) 


-1UDD1. 1 


u . u 


n 
u 


-ti, 1 


{■>) 


1 ACM 7 


A 7 A 


_i_n a ^of<- n 


E, 1 


(2) 


-10580.9 


80.2 


0; 


Ta, 1 


(3) 


-10505.6 


155.5 


±2.3112; 


T a , 2 


(3) 


-9251.5 


Aei 


±0.6634; 


E, 2 


(2) 


-9224.4 


Aei+27.1 


0; 


Tx, 2 


(3) 


-8161.7 


Ae 2 


±2.2447; 


T 2 , 3 


(3) 


-8154.2 


Aea+7.5 


±0.4503; 


Ai, 2 


(1) 


-8145.9 


Ae 2 +15.8 





E, 3 


(2) 


-8137.5 


Ae 2 +24.2 


0; 



by a contraction of the crystal lattice which stays cu- 
bic. We have associated this phase transition with the 
7 — > a transition of Ce. In real space the Pa3 ordering 
implies the appearance of four distinct sublattices of sim- 
ple cubic structure (see Fig. 3 of Ref. |l7]). We label these 
sublattices which contain the sites (0,0,0) (a/2)(0,l,l), 
(a/2)(l,0,l) and (a/2)(l,l,0) by {n p }, p = 1 - 4, respec- 
tively. In principle, one can proceed as in Ref. and 
derive an effective mean-field Hamiltonian. Here we will 
start from the crystal in real space and consider the fol- 
lowing four quadrupolar SAF's corresponding to the four 
sublattices of Pa3: 

5 {ni} (fi) = -^(^(O) ± S 2 (0) + S 3 (fi)), (4.13a) 

S {n2} (n) = -^(Sm - Sa(fi) ± 5 3 (n)), (4.13b) 

S {n3} (n) = ^(Si(fi) - S 2 (Q) - S 3 (fi)), (4.13c) 

S {ni} (n) = -^(Sm ± Sa(fi) - S 3 (fi)). (4.13d) 



Here we use the short notations Sh = S, 



(t=2,T 2g ,k=l-3) 



for 



real spherical harmonics Y 2 ls , Y 2 lc and Y 2 S which belong 
to a three dimensional irreducible representation T 2g of 
Oh- (These spherical harmonics are proportional to the 
Cartesian components yz, zx and xy for k = 1 — 3.) 

Below we consider the inter- site quadrupole interac- 
tions V QQ (n,n' ) which involve only the functions Q4.13a| - 
d). (There are alsa-SAF's with I = 4 and 6 allowed by 
the Pa3 symmetr yLD b ut those lead to weaker multipole 
interactions, Eq. ( 2.1 2| ) . ) We then rewrite Eq. (4.1) for a 
case when n 6 {n{\ and ft' G {n P >} (p' =2, 3, 4): 



{I fd WU\n'y QQ ^n%r fd )n'\J fd )n = -£■ 



{c{f.i}(W«) <KW»} { c {n p , } (w 3a> ) <KW> ) } 



(4.14) 
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where, as before, a, a', (3, /?' run over / and d, and where 
the same exclusion rules (if a = / then f3 = d and vice 
versa) hold between a and f3 as well as between a' and 
j3' . Here we have 



,aa a a 



dr r 2 / dr 1 r 



/2 



x TZl(r)TZl,(r')vA A {n,n'; r,r'), (4.15) 

with v\A(n,n'; r,r') where n — (0,0,0), ft' — 
(a/2)(0,l,l) and A = (I = 2, T 2g ,k = 1). The coeffi- 
cients Ci np \ are defined as 



2{n p }(Wa) = (i a \S {np] \j a ). 



(4.16) 



We introduce the quadrupolar density operators for the 
(fd) electron system on each sublattice: 

Paa(n P ) = E I 7 /**) c {n v }(i a j a ) S i0j/3 (J fd \, (4.17) 
i,J 

where n p £ {n P }, p = 1 — 4. Here again a = 0, = d or 
a = d, (i = f . In terms of quadrupolar density operators, 
the quadrupolar interaction operator between two (fd) 
systems at site n\ G {n\} and ri p > € {n p >} reads 



V(n 1 ,n p >) = 



Pi*(n x ) P^Mp')- (4-18) 



The mean-field potential at site n\ is obtained by sum- 
ming V(fii,n p ') over the twelve nearest neighbors n p > of 
ni on the fee lattice and by approximating the quadrupo- 
lar densities at these nearest neighbor sites by their ther- 
mal expectation values. The thermal expectation value 
of Pua(™p) does not depend (i.e. is the same) on any site 
of a given sublattice and from the equivalence of the four 
sublattices it follows that it is the same on all sites of the 
fee lattice. We then write 



(4.19) 



where the superscript e stand for thermal expectation. 
The mean-field potential is then given by 

U MF (ft l ) = -4^7 QaQ ' a 'pi«{nx)p%, (4.20a) 



or, explicitly, 



U 



MF i 



(n 1 ) = -(\ FF pf e + \^p% e )pt(n 1 ) 



DF <3,e^ 



(A^ P g< e + A^pr)P^i), (4.20b) 



\ FD J2- e \ „Q 



where we have defined X aa a a = A^ aa a a and used F 
and D for (act), with a = f or d, respectively. The 
values of \ aa a a have been calculated beforeila Here we 
use the values A FF =2241 K, X DF = \ FD =6489 K and 
X DD =18793 K, calculated for a-Ce. Including the intra- 
site part U\i n tra and the crystal field VcF(n\), Eq. (4.1), 
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FIG. 1. Calculated evolution of the order parameter am- 
plitudes pp' e and p^' e with temperature. 



which are also present in the ordered phase, we obtain 
the full mean-field Hamiltonian 



H 



MF i 



(n 1 )=U M "(n l )+Vc F (n{)+U\ i 



(4.21) 



Finally, the expectation values (the order parameter am- 
plitudes) p®' e and p% e , Eq. ( [4.19| ), are found by solving 
the following mean-field equations 



Q . e _ Tr{p Q F (n 1 ) eM-H^^/T]} 
Pf Tr{exp[-ff MF (n 1 )/T]} 

Q , e _ Tr{p%{Hx) eM-H MF (ft 1 )/T}} 



Pd 



Tr{exp[-# MF (ni)/T]} 



(4.22a) 
(4.22b) 



It is convenient to rewrite these equations in the basis 
\Kfd) = \kfkd) where H MF is diagonal, 



(4.23a) 



p Q D e = ^E c fn 1 }(W e " A ' /d/T ' ( 4 - 23b ) 



K f 



with 



(4.23c) 



Equations ( 4.20b )-( 1.23c ) are solved self-consistcntly. 
First, we introduce nonzero expectation values p®' e and 
p% e in the mean- field Hamiltonian H MF . After this 
we diagonalize H MF and calculate new values of p^' e 
and /)„ e at a given temperature T according to Eqs. 
( 4.22a ,b) . Then w e use these values to improve the mean- 
field Hamiltonian (4.21) and etc. until the input and out- 



put values of p®' e and pfj e converge. The results of the 
numerical calculation are shown in Fig. 1. The proce- 
dure outlined above converges very slowly in the vicinity 
of 100 K, i.e. at the phase transition point. We have 
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TABLE IV. Lowest levels of the energy spectrum of H MF 
at T = and magnetic moments M z for a-Ce. Numbers in 
parentheses stand for degeneracy; the site group is Se = Cg x i. 



r, n 




in K 


(ei - ei), in K 


Mz, in (iB 


E, 1 


(2) 


-10888.5 


0.0 


±2.0683 


A, 1 


(1) 


-10721.4 


167.1 





A, 2 


(1) 


-10699.7 


188.8 





E, 2 


(2) 


-10542.2 


346.3 


±1.0116 


E, 3 


(2) 


-10427.1 


461.4 


±0.4882 



found the transition temperature T± =97 K and the or- 
der parameter discontinuities Pp(Ti) = —0.06956 and 
p^ e {Ti) = —0.03875. (From symmetry considerations 
it follows that the phase transitio n is of first order .113) 
At T — the averages in Eq. ( 4.19| ) are taken over 
the ground state doublet and we obtain Pp £ (T = 0) = 
-0.15462, pg e (T = 0) = -0.08485. The lowest five lev- 
els of H MF for this case are given in Table IV. Notice, 
however, that unlike the cry stal field Vcf the mean field 
potential U MF , Eq. fl4.20bj ), and the Hamiltonian H MF 
depend implicitly on temperature T since the order pa- 
rameter amplitudes pp* e and p^ e change with tempera- 
ture, Fig. 1, and the energy splittings (ej — ei) decrease 
with increasing T up to T±. Above the phase transition 
point the spectrum transforms discontinuously to that 
of 7-Ce, Table III. The double degenerate states (the 
representations E in Table IV) are due to time-reversal 
symmetryES For a further discussion of the energy spec- 
trum of H MF we refer to Appendix A, where we study 
the magnetic moments of the two electrons. 



V. DISCUSSION AND CONCLUSIONS 

This work is an extension of our previous model of 
the 7 — a phase transition in Ce based on the idea of 
quadrupole ordering I13u3 In addition to the inter- site 
quadrupolar couplings here we consider the multipolar 
intra- site (direct Coulomb and exchange) interactions 
between one localized 4/ electron and one delocalized 
5d electron taken instantaneously at a same cerium site. 
The intra- site interactions are treated exactly in the 
adopted 4/5d model. In 7-Ce we have calculated and 
analyzed the crystal electric field excitations, Table III. 
In a-Ce the Pai quadrupolar ordering sets in and drives 
the 7 — a phase transition. The quadrupolar order has 
been s tudie d in t he mean-field approximation, Eqs from 
( 4.20b| ) to (4.23c). We have calculated the phase transi- 
tion temperature (T\ — 97 K) and the evolution of the 
order parameter amplitudes (p^ e and p§ e , see Fi g. 1) b y 
solving self-consistently the mearjrfijeld equations (4.20b)- 



( 4.23c ). We have shown beforeOES that quadrupolar or- 
dering in the Pa3 structure leads to a uniform lattice 
contraction conserving the cubic symmetry of the lattice. 



The change of the energy spectrum at the transition 
implies a change of the magnetic susceptibility. Indeed, 
we have found that the calculated magnetic moments 
are different in 7— and a— Ce. Moreover, in the ground 
state of a— Ce the magnetic moment of the 4/ electron is 
bound to the magnetic moment of the 5c? electron (a qual- 
itative origin of this correlation is given in Appendix B). 
The lowest magnetic excited state (E, 2 in Table IV) 
is separated from the ground state by an energy gap 
Ae ~ 350 K which is much larger than a typical crys- 
tal field excitation in the 7 phase, Table III. However, 
here our treatment is incomplete. Although the present 
model carefully takes into account the intra- site interac- 
tions it does not describe properly the metallic bonding 
in Ce. 

The question of correspondence between our approach 
and electron band structure calculations deserves a spe- 
cial attention. As we have discussed in Sec. II, in the 
"muffin-tin" approximation a localized 4/ electron expe- 
riences only a field of spherical symmetry and occupies 
a 14-fold degenerate level. The localized states of the 4/ 
electron then are uncorrelated with the states of conduc- 
tion electrons, because the spherical component of the 4/ 
density is independent of its spin (s z ) and orbital (rrif) 
projections. In our study we show that this simple pic- 
ture is not correct and there exist strong local correlations 
between localized 4/ and delocalized 5c? electrons omitted 
in a conventional band structure calculation. These cor- 
relations arise due to the Coulomb on-site repulsion and 
reflect the electronic term structure of atomic cerium, 
Sec. III. We show that the excitations of the 4/ electron 
are combined with those of 5c? electron in a single spec- 
trum, which is sensitive to crystal site symmetry because 
of intersite interactions, Sec. IV. 

In principle, band structure calculations with the full 
potential (FPX. extension (so-called FP-linear muffin- 
tin (FLMXO)c3 and FP-linear augmented plane wave 
(FLAPW)E3 methods) are capable of dealing with non- 
spherical contributions of density and potential. Pro- 
vided that the site symmetry is introduced explicitly, 
calculations with the full potential option can describe 
some, but not all structural properties associated with 
the Fm3m — ► PaS transformation. The reason is that 
the band structure calculations are based on the single- 
determinant Hartree-Fock method. In our treatment 
each local two-electron basis function, Eqs. fl2.2| )-(2.6), 
corresponds to a Slater determinant (with the permuta- 
tion property (2.7)). The solutions are expressed as lin- 
ear combinations of all these functions (determinants). 
As such, our method corresponds to a many determinant 
treatment or configurational interaction (CI). Clearly, it 
is not the case with the band structure approach. This 
explains why the conventional band structure treatment 
is missing some intra- and inter- site correlations which 
are taken into account in our approach. 

The other drawback of the FP electronic band struc- 
ture calculations is connected with the local density ap- 
proximation (LDA). In the LDA the exchange potential 
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is a function of density and has the full (unit) symme- 
try of the crystal. This implies that the density and the 
LDA exchange potential are invariant under inversion in 
both phases. Expanding the exchange potential in terms 
of spherical harmonics on a cerium site we find that only 
harmonics with even I are allowed by the inversion sym- 
metr y, i.e., I = 2, 4 etc. However, as we observe from 
Eq. ( 3.10b ) the contributions with odd values of I are 
relevant for the exchange between 4/ and 5d electrons, 
namely I — 1, 3 and 5. Thus, these terms are washed 
away by the LDA treatment. 

On the other hand, the present method should be 
extended to include the metallic bonding in Ce explic- 
itly. In our opinion, a combination of the local corre- 
lations with the band structure approach constitutes a 
challenge for further studies. We are presently working 
on the problem and hope to achieve this by employing 
the valence bond (VB) or Heitler-London approach. The 
VB method is more difficult to implement, but usually 
it gives a better description of .djicmical bonding than 
the method of molecular orbitals£3a The alternative ap- 
proach is a merger with one of the existing band structure 
methods and introducing in some way a CI treatment. 

Finally, we would like to mention again that we are 
aware of the fact that our approach is not complete, but 
it certainly underlines the importance of the structural 
factors and the local correlations for this long-standing 
problem. We have shown that the local electronic inter- 
actions can trigger the 7 — > a phase transition in Ce, 
but other aspects of the problem (in particular, chem- 
ical bondinft_etc.) should also be taken into account. 
As beforeJiMj we suggest synchrotron radiation experi- 
ments in order to check the appearance of weak super- 
structure reflections in the Pa3 structure (a-Ce) and to 
study diffuse scattering in 7— and a— Ce. 

The present model generalizes our initial approach, 
Rcf. [l7], for a many electron case. The method can be 
applied (as it has been done, for TmTe in Ref. ^6|) to 
study quadrupole ordetingsL^ m, lanthanides, and theit 
compounds (DxB 2 C 2 ,eZl DyB 6 ,E2 UCu 2 Sn,N PrPb 3 ,El 
YbAsj!II YbSbE3) where a few / electrons at each site 
are involved in the process of ordering. 
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APPENDIX A: 

In order to calculate the effective magnetic moments 
and the Lande g-factors we study the polarization of elec- 
tronic states in a small external magnetic field H . In such 
case we add to the Hamiltonian TL (that is U\i n t ra , Eq. 



( |3.16 ), for the atomic case, U \ intra + Vc e f, Eq . ( 4.11 ), for 



7-Ce, and U\ intra + VcEF 
a magnetic term 



U MF , Eq. fl4.2l|) , for a-Ce) 



Here 



with 



M z 



V mag = -M z ■ H. 



(M z (f) + M z {d)l 



M(f)=ML(f) + 2S(f)), 
M(d) =fi B (L(d) + 2S(d)), 



(Ala) 

(Alb) 

(A2a) 
(A2b) 



fiB being the Bohr magneton. The matrix elements of 



Vmag read 



(Ifd\V mag \ Jf d ) 



((if\M z (f)\j f )S Uji 
(i d \M z {d)\j d )S ifjf )-H, 



(A3) 



where (if\M z (f)\jf) and (i d \AA z (d)\j d ) are one-particle 
matrix elements which can be easily computed. If we di- 
agonalize the matrix of {H. + V ma g), then the degeneracies 
of the energy terms are lifted and the magnetic moment 
of each sublevel v is given by 



M z (e v ) = (e v \M z \e u ), 



(A4) 



where {e„} stands for the energy levels {Efd} or {ei} 
given in Tables II or III and IV, for the atomic case (intra- 
site), the 7-phase and the a-phase, respectively. The 
results for A4 Z for the solid state phases are quoted in 
Tables III and IV. The Lande g-factors are obtained from 



M z (e v ,j) = fi B g(£v)j, 



(A5) 



where j is the z-projection of a total angular moment: 
j = —J.—J+ 1, ...,+J. The calculated (/-factors for 
the atomic case, Table II, are close to the experimental 
values £.3 

In case of 7-Ce and a-Ce the z-axis is the [001] axis of 
the cubic crystal. For 7-Ce only triply degenerate levels 
of T\ or T 2 symmetry have nonzero magnetic moments 
which are — A4(t), and +M(t), where M{t) > and 
t = (Ti,/x) or t = (T 2 , fi), Table III. For a-Ce only double 
degenerate levels of E symmetry have nonzero magnetic 
moments which are — M(e) and +M(e), where M(e) > 
and e = (E,/j,), Table IV. The double degenerate states 
are due to time-reversal symmetry£a It is worth noting 
that in the ground state (E,l) the magnetic moment of 
the localized 4/ electron is not independent. It is at- 
tached to the magnetic moment of the 5d conduction 
electron. The change of sign of the magnetic moment of 
the 5d electron under the time reversal symmetry now 
requires the concomitant change of the sign of the mag- 
netic moment of the localized 4/ electron. We discuss 
the origin of this correlation in Appendix B. Since the 
5d electron belongs to the conduction band, the ground 
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FIG. 2. The transformation {x,y,z} — > {x',y',z'} (the 
Euler angles are a = 7r/4, /3 = arccos (l/v3), 7 = for 
passive and (7, /3 , a) for active rotations), and the function 
5 {ni} , Eq. fl4.13aP. 



state (E,l) gives rise to a Pauli paramagnetism. The 
other levels given in Table IV represent excited states of 
a-Ce. However, we observe that the first (representation 
A,l of Table IV) and the second state (A, 2) are nonmag- 
netic. The first magnetic excited state (E,2) is separated 
from the ground state by an energy gap of Ae ~ 350 K. 
This observation could explain the absence of the Curie- 
Weiss contribution to the magnetic susceptibility \ a at 
temperatures T < T\. 



APPENDIX B: 



We consider the function <S/ ni }(fi), Eq. (4.13a), in a 
new coordinate system of axes {x',y',z'} where the z'- 
axis corresponds to the cubic [lll]-axis, Fig. 2. Then 



5 {ni} (0) 



Y, 



m=0 
1=2 



(Bl) 



where il' stands for the polar angles (Q',(/>') in the 
new system of axes. The coefficients cf n ^-,(ifjf) and 

c® ni y(idjd)i Eq. (4-16), become diagonal in the new basis, 



C{ ni }(i a ja) = (i a \Y2\i a ) Si 



(B2) 



that facilitates the c alculat ion of the mean-field interac- 
tion U MF (n 1 ), Eq. ( |4.20h| ). In this Appendix we show 
that JJ is minimized for the following two functions 



y 



Mi =5 
L=5 



L=5 



Yi(n' f )-Y*(n' d ), 
= Y 3 - 3 (n')-Y 2 - 2 (n' d ), 



(B3a) 
(B3b) 



where fl'j and Sl' d refer to the polar coordinates of 4/ and 
5d electron, respectively. 
Indeed, there are 11 



functions 



y A L % (m l = 

—5,-4, ...,5), which form a basis of the two-particle ir- 
reducible representation H of the group SO(3) of three 
dimensional rotations. (In Eq. (B3) we quote only two 
functions, the others can be obtained from the Table 



4 3 of Ref. E^.) In the LS (Russel-Saunders) coupling 
the triplet 3 ~ff lies lower in energy than the singlet 1 H, 
Table I. If now we include the diagonal mean-field cou- 
pling U MF ', we obtain that the six states with the orbital 
components and 3^~ 5 and with the spin components 
Ms = —1, 0, 1 (originating from the 3 H triplet) go down 
in energy. These six states are further split by the spin- 
orbit coupling in three magnetic doublets of E symmetry 
of the site group 5 6 = C^xi. Finally, in the full treatment 
(Sect. 4.2) the middle doublet is split by the crystal field 
in two nonmagnetic components of A symmetry. This 
explains qualitatively the origin of the four lowest levels 
of the Pa3 structure of a-Ce, Table IV. The real situa- 
tion is more complex since there is a mixing of different 
configurations in atomic cerium. Nevertheless, a consid- 
erable admixture of 3 H configuration (29%) is found in 
the ground state of atomic ceriumEj and it is very likely 
to occur in 7-Ce. 

In a-Ce orbital functions of 5d electrons on neighboring 
sites overlap giving rise to band structure effects. Here 
we consider these effects only as a perturbation to the 
ground state E, 1, Table IV. Then in tight- binding ap- 
proximation a band state with the wave vector A; is a 
mixture of the d— functions V 2 2 (f2^) and (Q^). How- 
ever, as we have seen earlier, these electronic 5c? states 
are bound to the two 4/ electron states 1^ 3 (0^) and 

y-3 



(fl'f) and form the two-electron functions J^ffg 



and y^=f & , Eqs. (|B3a| ,b). Therefore, the resulting 
doubly degenerate states are in fact two-electron band 
states with the energy £+(k) = £-(k) (the sign + or 
— here stands for the two-electron time reversed band 
states). In the magnetic field these two bands become 
split, £+{k) 7^ £_(fc), and this leads to a temperature 
independent Pauli paramagnetism. 
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